Effect of uniaxial strain on the structural and magnetic phase transitions in BaFe2As2 
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We report neutron scattering experiments probing the influence of uniaxial strain on both the 
magnetic and structural order parameters in the parent iron pnictide compound, BaFe2As2. Our 
data show that modest strain flelds along the in-plane orthorhombic b-axis can affect significant 
changes in phase behavior simultaneous to the removal of structural twinning effects. As a result, 
we demonstrate in BaFe2As2 samples detwinned via uniaxial strain that the in-plane C4 symmetry 
is broken by both the structural lattice distortion and long-range spin ordering at temperatures far 
above the nominal (strain- free) , phase transition temperatures. Surprising changes in the magnetic 
order parameter of this system under relatively small strain fields also suggest the inherent presence 
of magnetic domains fiuctuating above the strain-free ordering temperature in this material. 



PACS numbers: 74.70.Xa, 74.62. Fj, 75.50.Ee, 75.40.Cx 

A number of recent experiments have exploited uni- 
axial strain as a means of structurally detwinning bi- 
layer iron pnictide compounds in order to explore the 
anisotropics inherent in the electronic ground states of 
these materials [IHl] ■ A key result of many of these stud- 
ies in both undoped and doped bilayer iron pnictides has 
been the report of dramatic anisotropics in the in-plane 
charge transport that surprisingly persist to tempera- 
tures far above the strain-free tetragonal to orthorhombic 
phase transition in these materials (Tg) [I]. This in-plane 
transport anisotropy develops between the orthorhombic 
a- and 5-axes and is indicative of broken C4 symmetry in 
the tetragonal phase induced via the applied, detwinning, 
strain field. The surprising evolution of this anisotropy as 
a function of doping has led to a novel picture involving 
an additional energy scale in the iron pnictide phase di- 
agram possessing C2 symmetry associated with nematic 
order [6]. 

Within the detwinned parent AFe2As2 (A = Ca, Ba, 
Sr, ...) (A-122) compounds, the only systems to date that 
exhibit this anisotropic in-plane behavior above T5 are 
BaFe2As2 [1] and EuFe2As2 [5], both of which possess 
more continuous structural and magnetic phase transi- 
tions relative to other 122 variants ^(8^. This peculiar- 
ity suggests that the nominal (strain-free) phase behavior 
in these systems — and in particular BaFe2As2 — may be 
particularly susceptible to the the influence of mechan- 
ical strain relative to other 122 counterparts whose in- 
trinsically stronger first order phase behaviors preclude 
anisotropy resolvable above their strain- free T5 [21 H]. 

In Ba-122, where the phase behavior is seemingly the 
most continuous, an anomalous lattice softening has been 
observed in resonant ultrasound experiments in proxim- 



ity of T5 [iTl HHj thus suggesting a large susceptibility 
to the influence of applied strain fields. This fact com- 
bined with the growing number of studies utilizing uni- 
axial strain as a means of detwinning Ba-122 and re- 
lated 122-compounds [S] suggests that it is particularly 
important to directly explore the influence of the result- 
ing strain flelds on both its magnetic and structural phase 
transitions. 

In this letter, we report neutron scattering measure- 
ments probing the effects of uniaxial, compressive, strain 
on both the magnetic and structural phase transitions 
within an iron pnictide parent compound, BaFe2As2. 
Neutrons allow for the bulk magnetic and structural re- 
sponse of the sample to be surveyed simultaneously as 
the planar strain held is enhanced, and phase behavior 
in the system can therefore be directly observed as the 
strain fields necessary to detwin a given crystal are ap- 
plied. Similar to other iron pnictide parent materials |101 - 
IT2] Ba-122 develops an ordered antiferromagnetic phase 
(Taf) with Qaf = (7r,0,7r) below T5; however, within 
Ba-122 and other parent 122 systems these two transi- 
tions develop nearly simultaneously [111 [15] . The pres- 
ence of strong spin lattice coupling in these 122 parent 
pnictides [HKISldS] along with the potential influence of a 
tricritical point [TS] in their phase diagrams render sharp 
phase transition onsets corresponding to either strongly 
first order [16| or nearly second order 8] phase transitions 
whose nature varies with the A-site ion. 

Our data show that as compressive uniaxial strain 
fields are applied along the orthorhombic 6-axis, the nom- 
inally sharp structural phase transition shifts upward and 
broadens in temperature accompanied by a surprising 
parallel increase in the onset temperature of long-range 
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FIG. 1. (a) Picture of detwinning mount used for neutron 
scattering experiments, (b) Reference data (from Ref. 8) re- 
fit with a single Gaussian fit quantifying the structural phase 
transition in strain-free Ba-122. (c) Single Gaussian fits to 
radial scans through the (2, 0, 0) reflection showing a broad- 
ening below Tg using data from Ref. 8. 



antiferromagnetic order. Our results suggest that the 
strong in-plane anisotropies observed above the strain- 
free T5 in detwinned Ba-122 can be primarily accounted 
for by a direct shift and renormalization of the magnetic 
and structural order parameters under applied strain. 

BaFe2As2 single crystals were grown via standard, self- 
flux techniques [19]. Once grown, the samples were cut 
along the orthorhombic a- or 6-axes and mounted within 
a small pressure clamp (Fig. 1 (a)). A belleville washer 
was used to monitor applied force across the sample. 
Samples were aligned within the [H 0, L] scattering plane 
using the Fmmm orthorhombic unit cell. Neutron exper- 
iments were performed on the HB-lA triple axis spec- 
trometer at the High Flux Isotope Reactor, Oak Ridge 
National Lab and on the N5 triple-axis spectrometer at 
the Canadian Neutron Beam Centre, Chalk River Lab- 
oratories. Experiments on N5 were performed with a 
pyrolitic graphite (PG) monochomator {Ei — 14.5 meV) 
and analyzer with a PG filter placed after the sample and 
coUimations of 30' - 60' - sample - 33' - 144' . The HB- 1 A 
setup consisted of a double-bounce PG monochromator 
(Ei — 14.64 meV), PG analyzer with a PG filter before 
the sample collimations of 48' — 48' — sample — 40' — 68'. 

Data probing the evolution of the structural order 
parameter as pressure is increased along the [0, 1, 0]- 
orthorhombic axis are plotted in Figure 2. The metric 
by which the structural order parameter can be resolved 
necessarily evolves as the pressure is increased and the 
crystal is progressively detwinned. Under zero strain, 
the formation of twins below T5 in principle allows both 
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FIG. 2. Peak widths and lattice parameters determined from 
Gaussian fits to the Q = (2, 0, 0) nuclear reflection plotted 
in panels (a) and (b) respectively. Fit parameters show the 
evolution of the phase behavior as pressure along the 6-axis 
is increased incrementally in the same crystal. Insets show 
global picture of the phase transition. Radial scans through 
Q = (2, 0, 0) both above and below Tg with (c) MPa and (d) 
0.7 MPa of uniaxial pressure applied. Horizontal bars show 
the full width at half maximum (FWHM) of the Gaussian fits. 



the a- and 6-axes to be resolved within the same scatter- 
ing plane and the structural order parameter S = to 
be determined; however the momentum resolution of our 
neutron experiments only resolves a broadening of the 
peak line shape below T5 as the orthorhombic a- and 
&-axes distort away from each other in twinned samples. 
Due to this and in order to obtain an unbiased fit to 
the structural peaks as various strain fields are applied, 
we instead fit all peaks as single Gaussians whose width 
changes as the peak splits and twin domains develop. 

This single Gaussian treatment of the structural peaks 
is appropriate in the case of a detwinned crystal; how- 
ever it is only an approximation for two closely spaced 
Gaussian peaks in fits to data from twinned crystals. To 
verify that this single peak approximation is adequate 
in the twinned limit, we remodeled previously reported 
data (using similar collimations and the same spectrom- 
eter) using only a single Gaussian fit rather than two 
symmetrically displaced Gaussian peaks and plotted the 
results in Figs. 1 (b) and (c). Refitting the data gives 
an identical onset temperature of Taf = 136 K to that 
previously reported fSJ demonstrating the validity of this 
approximation for determining transition temperatures. 

Turning now to the current BaFe2As2 data, the peak 
widths and resulting lattice parameters determined by 
the Gaussian peak fits to the (2, 0, 0) nuclear Bragg peak 
are plotted as a function of temperature for three strain 
fields in Figs. 2 (a) and (b) respectively. The (2, 0, 0) 
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peak width broadens at Taf = 135.7 K in the strain-free 
case consistent with the expected Tg [8] ; however as pres- 
sure along the [0, 1, 0] axis is increased progressively to 
0.7 MPa the broadening observed at Tg vanishes. This 
is exphcitly shown in Figs. 2 (c) and (d) where the (2, 0, 
0) peak is shown above and below T5 in both MPa and 
0.7 MPa respectively. Within the resolution of our mea- 
surements, there is no resolvable width change in peak 
line shape once 0.7 MPa is applied suggesting that the 
crystal has been substantially detwinned. 

In order to resolve the onset of the structural phase 
transition as the crystal is progressively detwinned, the 
longer in-plane lattice parameter a is plotted as a func- 
tion of temperature in Fig. 2 (b) for the three applied 
pressures. Here, the shift in the average in-plane lattice 
parameter for the strain-free crystal simply reflects the 
asymmetric splitting between a- and &-axes as the struc- 
tural phase transition sets in; however under 0.7 MPa 
of pressure, it is immediately apparent that the onset 
of T5 has shifted upward in temperature to Tg « 147 
K. Fig. 2 (b) explicitly demonstrates that as pressure is 
increased, the onset temperature for T5 shifts systemati- 
cally upward as strain fields approaching the detwinning 
threshold are approached. 

Now turning to the response of spin ordering to the 
applied uniaxial strain field, the intensity of the Q — 
(1,0,3) magnetic peak is plotted in Fig. 3 (a). Similar 
to the response of the structural phase transition, the on- 
set of magnetic order systematically increases as strain 
is increased across the sample. Specifically, an enhanced 
tail of magnetic scattering appears that deviates from the 
previously modeled power law behavior [8j. Our previ- 
ous measurements in strain-free BaFe2As2 have already 
demonstrated that this tail of scattering is present across 
a similar range of temperatures (albeit at much lower in- 
tensities [S]) and suggest that the application of strain 
enhances the volume fraction of the ordered moment for- 
mation across this higher temperature scale. The spin 
order that develops in this tail of scattering shown in Fig. 
3 (c) is long-range with a minimum spin-spin correlation 
length at 142 K of ^ = 208 ± 19A— identical within error 
to data collected at 3K in this sample. The diameter of 
the correlation length was found via ^ = 2v/2ln(2)(w)-i 
where w is the fit Gaussian peak width in A~^. 

In Fig. 3 (a), at T = 3 K the relative scattering in- 
tensity from magnetic moments contributing to the (1, 0, 
3) reflection (determined by normalizing the (1, 0, 3) to 
the (0, 0, 4) nuclear peak intensity) increases upon the 
application of increased strain. This is simply due to the 
progressive reduction of the volume fraction of the sam- 
ple with magnetic domains whose moments are oriented 
out of the scattering plane. We note however that if we 
assume no twinning under 0.7 MPa the ordered moment 
at 3 K is only 0.87/zb relative to IMfis measured in the 
twinned, zero strain, state — suggesting that the sample 
remains partially twinned (sa 43%) at 0.7 MPa pressure. 
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FIG. 3. (a) Intensity (/ oc M^) at the Q = (1, 0, 3) AF peak 
position normalized by the (0, 0, 4) nuclear peak plotted as a 
function of temperature for different pressures applied along 
the 6-axis. Inset shows magnetic phase transition over the 
entire temperature range. Dashed lines are power law fits 
as described in the text, (b) Pressure dependence of Taf, 
Ts, and TAF.onset- Open symbols denote measurements on a 
seperate crystal at an intermediate pressure, (c) Scan Q=(l, 
0, 3) at T= 142 K with nonmagnetic background subtracted. 
Horizontal bar indicates the FWHM of the (1, 0, 3) reflection 
measured at 3 K. (d) HB-IA data measuring the magnetic 
peak intensity (circle symbols) at Q=(l, 0, 3) under an un- 
known P> 1 MPa. Squares show the a-axis lattice parameter 
measured under the same strain at the Q=(2, 0, 0) nuclear 
position. Triangles show the same a-axis measurement with 
the strain released. Solid lines are mean field fits. 



The tail of scattering above Taf in Fig. 3 (c) at 
present does not seem to originate from a simple picture 
of critical fluctuations above a second order phase tran- 
sition given the lack of a resolvable divergence in the 3D 
spin-spin correlation length; rather the enhanced mag- 
netic scattering seems to stabilize over large domains 
within the sample as strain is increased. Nevertheless, 
critical scattering can not be entirely precluded so we 
have instead identified three energy scales: (1) an AF or- 
dering temperature Taf determined with a simple power 
law = A{1 — T/TafY^ using the range of temper- 
atures 120 K< T < 136 K, (2) an onset temperature, 
TAF,onset j whcrc three dimensional magnetic scattering is 
first resolvable, and (3) the onset of the structural phase 
transition T5 . These three energy scales are plotted as a 
function of increasing pressure in Fig. 3 (b). 

The structural phase transition shifts upward approx- 
imately 11 K upon applying modest uniaxial pressure of 
only « 0.7 MPa. Increases in both Taf and T5 manifest 
with strain fields below the detwinning threshold; how- 
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ever, upon applying enough strain to substantially de- 
twin the sample, both TAF,onset and Tg increase dramat- 
ically. Our results are therefore germane to studies prob- 
ing intrinsic anisotropics within BaFc2As2 under uniaxial 
strain (6j. Rather than simply allowing the resolution of 
intrinsically anisotropic electronic fluctuations that break 
the C4 symmetry of the tetragonal phase, strain fields 
along the [0, 1, 0]-axis of Ba-122 generate a shift in the 
nominal, strain-free, phase behavior. From our current 
data however, it is not clear which shifted phase primar- 
ily drives the observed transport anisotropics in strained 
Ba-122. The surprising rise in the onset of AF order with 
the application of relatively small strain fields also sug- 
gests the influence of two possible effects in Ba-122: (1) 
the presence of magnetic domains fluctuating far above 
the strain-free Taf |20j that are eventually pinned by the 
strain fleld and stabilize at higher temperatures as T5 is 
simultaneously shifted upward and (2) the influence of 
strong spin-lattice coupling in this system resulting in an 
enhancement of the magnetic order via direct coupling 
to the applied strain field. 

In order to test this notion, we applied pressure on 
a separate sample beyond the tolerance of the belleville 
washer in the sample mount and generated a 6-axis pres- 
sure P> 1 MPa, and the resulting magnetic order pa- 
rameter and a(T) are plotted in Fig. 3 (d). The onset 
of Tfi and Taf appear to decouple under this pressure 
with Ts « 155 K and T^f ~ 145 K (determined from 
a mean field fit = ^(1 - T/Taf)), and the mag- 
netic data show that for sufficiently high strain fields the 
enhanced tail of the magnetic order parameter fills in al- 
most entirely. The magnetic order parameter seemingly 
approaches a truly broadened magnetic order parameter 
with no clear distinction between a broadened magnetic 
tail and the residual bulk onset of AF order and does not 
simply continue to shift upward in temperature. 

This data when combined with the known presence of 
a muted tail of magnetic scattering in strain-free Ba-122 
across an identical temperature range [8. suggest that the 
stabilization of fluctuating magnetic domains constitutes 
the primary origin of the increase in T AF,onset- Rather 
than a continual increase in the temperature scale of the 
magnetic tail via spin lattice coupling, instead the mag- 
netic tail flUs in as strain is increased and T5 is driven 
further upward allowing the intrinsic phase behavior of 
the magnetic order to manifest itself when T5 is driven 
above « 147 K. As a further check, we also measured 
the magnetic order parameter in a separate Ba-122 crys- 
tal annealed for several days in vacuum (where transport 
anisotropy is substantially diminished|22|): however this 
failed to produce an enhanced tail or comparably changed 
magnetic order parameter — demonstrating that uniaxial 
strain does not simply serve to mitigate the influence of 
defects within the Ba-122 matrix. 

In conclusion, we have investigated the influence of uni- 
axial strain on the phase behavior of the parent bilayer 



iron pnictide system BaFe2As2. Contrary to prior as- 
sumptions [6], we observe a progressive increase in the 
onset temperatures of both Tg and Taf with increas- 
ing uniaxial strain along the in-plane b-axis. Once the 
crystal has been substantially detwinned, both the onset 
of the magnetic and structural phase transition temper- 
atures have shifted upward dramatically and their phase 
behavior has broadened into a smeared, second-order-like 
transition. Our results suggest that the in-plane charge 
transport anisotropy reported above the strain- free Tg 
in mechanically detwinned BaFe2As2 [H H] likely stems 
from either or both of the renormalized magnetic and 
structural phase transitions under enhanced strain flelds. 

Note added: Upon preparing this manuscript for sub- 
mission, a related work was posted probing the struc- 
tural phase transition's response to uniaxial strain in 
BaFe9As9 f21j. Their result reporting an upward shift in 
T5 upon the application of strain is consistent with the 
results presented here. 
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